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Phonon statistics in an aconstical resonator coupled to a pumped two-level emitter 

Victor Cebar0 and Mihai A. Macovelll 

Institute of Applied Physics, Academy of Sciences of Moldova, Academiei str. 5, MD-2028 Chiyindu, Moldova 

(Dated: February 5, 2015) 

The concept of an acoustical analog of the optical laser has been developed recently both in 
theoretical as well as experimental works. We discuss here a model of a coherent phonon generator 
with a direct signature of quantum properties of the sound vibrations. The considered setup is made 
of a laser driven quantum dot (QD) embedded in an acoustical nanocavity. The system’s dynamics 
is solved for a single phonon mode in the steady-state and in the strong QD-phonon coupling regime 
beyond the secular approximation. We demonstrate that the phonon statistics exhibits quantum 
features, i.e. sub-Poissonian phonon statistics. 

PACS numbers: 42.50.Ct, 42.50.Lc, 78.67.Hc, 43.35.Gk 


I. INTRODUCTION 


Sub-Poissonian statistics of vibrational states is a pure 
quantum propriety of an oscillating system. The domain 
of this statistics starts at the limit of a classical coherent 
state having a Poissonian distributed quanta and may 
end up with a pure Fock state at the other limit. Stud¬ 
ies on the quanta statistics had already revealed many 
pure quantum features for different physical systems and 
remarkable results were achieved in a large spectrum of 
jhotonic quantum electrodynamics’s (QED) applications 
11 as well as in Bose-Einstein condensate’s (BEC) physics 
2,1^. More recently, successful experiments in cooling 
and detection of mechanical systems in the near ground 
state domain enhanced a particular interest in the 
research of similar features in the acoustical domain, due 
to bosonic quantification of the sound’s vibrations. 


Earlier experiments in laser generation of coherent 
phonons in different bulk materials were succeeded 
by new optomechanical and electromechanical setups 
in phonon QED, achieving important experimental re¬ 
sults in the acoustical analog of the optical laser by us- 
ing piezoelectrically excited electromechanical resonators 
and laser driven compound microcavities or 
trapped ions [T^ . In the meantime, theoretical mod¬ 
els propose improvements in the background theory of 
the experiments like the PT-symmetry approach [l3l| and 
two cavity optomechanics [3 j as well as new possible se¬ 
tups using vibrating membranes [T3 . [^I, quantum dots 
embedded in semiconductor lattices Il6l-fl3 and BECs 
under the action of magnetic cantilever [l9j. Moreover, 
quantum features like sub-Poissonian distributed phonon 
fields have been already predicted in optomechanical se¬ 
tups based on vibra ting mirrors and in single¬ 

electron transistor s 12211 as well as phonon antibunch¬ 
ing [T^ , squeezin g |20l l23| and negative Wigner function 
of phonon states [24| . In addition to increasing perfor¬ 
mances of the optomechanical devices [1^ , reports on the 
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art-of-the-state of acoustical cavities [26l - l^ have shown 
good phonon trapping in bulk materials with high cavity 
quality factors, thus leading to the concept of acoustical 
analog of photon cavity-QEDs (cQED). 

In this article, we study the model of a coherent phonon 
generator in a setup consisting of a qubit embedded in 
an acoustical cavity involving strong QD-phonon-cavity 
coupling regime. The qubit, i.e. a two-level quantum dot, 
is driven by an intense laser field and acts as a phonon 
source. Under action of laser light the electron jumps 
from the QD’s valence band to the conductance band 
leaving a hole in the valence band. The created exci- 
ton (electron-hole) represents the QD’s excited state and 
interacts with the acoustical vibrations, thus, creating 
or annihilating phonons in the cavity. We demonstrate 
that in analogy with recent experiments in photon cQED 
[2^ the strong qubit-resonator couplings may reveal ad¬ 
ditional quantum phenomena in the phonon cQED’s do¬ 
main. Particularly, we show that in this regime the gen¬ 
erated steady-state phonon held obeys sub-Poissonian 
phonon statistics. 

This paper is organized as follows. In Sec. II a detailed 
description of the model is given, i.e., after an introduc¬ 
tion of the used nomenclature we focus on simplifying the 
system’s Hamiltonian in order to arrive at an easy solv¬ 
able master equation for the reduced density operator of 
the QD-phonon system. In Sec. Ill, we present a general 
overview over the model’s results and we further discuss 
the important aspects of the study. The Summary is 
given in Sec. IV. 


II. THEORETICAL FRAMEWORK 

A two-level laser-pumped semiconductor quantum dot 
is embedded in an acoustical nanocavity (see also [l^ . 
E3)- The QD’s transition frequency between its ground 
state \g) and the excited state |e) is denoted by ujqd- The 
excited QD may spontaneously emit a photon with 7 
being the corresponding decay rate, see Fig. [TJ For a 
more realistic case, we introduce the dephasing losses rate 
through 7c. The single mode cavity phonons of frequency 
ojph are described by the anihilation (b) and the creation 
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FIG. 1: (color online) The schematic of the investigated 
model: A two-level QD is fixed in a mnltilayered strncture 
forming the acoustical nanocavity. The QD is pumped near 
resonance with a coherent laser sonrce of frequency ljl- The 
emitter may spontaneously emit a photon at a decay rate 7 
while the cavity phonon damping rate is denoted by k (see 
also Ref. Hi). 


[b^] operators, respectively. The system’s Hamiltonian 
is: 

H = nujgdS, + nujphb'<b + hn{S+e-^‘^^* + 

+ hgS+S-{b^+ b), (1) 


follows: 

H — Hslow 

Hfast = hg{sm^ OR _-h cos^ 0i?++){6'*'e“p'‘* -I- H.c.} 

_ ^^ sm ( 20 ) ^ 4 ) 

Instead of making the usual secular approximation [sO, 
lUI we do keep the fast rotating terms in the QD-phonon- 
cavity interaction Hamiltonian. Their main contribution 
is evaluated as [H, [s^ : 

^fast ~ ~J^Hfast{t) j dt Hfast{t) 

= Ho — hARz + hpb^bRz, (5) 

with 

^ ^ ^ / cos ( 20 ) _ sin^ ( 20 ) \ 

2 \ oJph 4:{ujph + 2Q) J 

and 


where the QD’s operators are defined as: 5+ = \e){g\, 
S~ = \g){e\ andS'^ = \ (|e)(e| - \g){g\) obeying the stan¬ 
dard commutation relations for SU(2) algebra. The first 
two terms correspond, respectively, to the unperturbed 
QD and to the free single mode phonon Hamiltonians. 
The third term corresponds to the QD-laser interaction 
within rotating wave and dipole approximations whereas 
Wi is the laser frequency. The last term describes the 
QD-phonon-cavity interaction with g being the coupling 
strength constant. 

In what follows, we describe the Hamiltonian in a 
frame rotating at the laser frequency ujl and apply the 
dressed-state transformation: 

|-|-) = sin0|(/)-I-cos0|e), 

|-) = cos 0 |g) - sin 0 |e), ( 2 ) 

where 20 = arctan(2Q/A) while A = ujqd — ujl is the 
detuning of the laser from the QD’s transition frequency. 
The dressed-state system’s Hamiltonian becomes then; 

H = hflRz + hujphb^b + hg{b^ + b){sm^ 6 R _ 

-b cos2 0i?++-^(i?+-bi?-)}, (3) 


= 9 


sin^ (20) 
4(wp/j -b 2 fi) 


Here, Hq is a constant and can be dropped as it does 
not contribute to the system’s dynamics. Notice that 
the coupling regimes are related not only to the QD- 
phonon-cavity coupling constant g but also to the con¬ 
tribution coming from fast rotating terms, i.e. A and /?, 
proportional to g^. Consequently, for weak QD-phonon 
couplings strengths g the fast rotating terms contribu¬ 
tion can be neglected, i.e., {A,/3} = 0. For strong QD- 
phonon coupling regimes the secular approximation is no 
longer justified and the contribution of plays a role 

which will be considered. Thus, the final Hamiltonian 
H = Hslow + Hfait 


H = h{ujph — 2Q)b^b — hARz -b h(3b^bRz 

- hg^-^^^{b^R- +R+b). ( 6 ) 


To solve the QD-phonon system’s dynamics one uses 
the density matrix formalism for the reduced density op¬ 
erator p\ 


P 


^ p] ~b diqdP ~b Ap/iP, 


(7) 


where D = -|_ (A/2)2. The new QD’s operators are: 

R^ — |+)(—I, R~ = |—)(+|, R++ = |+){+|, R — = 

|—)(—I and Rz = i?++ —i?__ satisfying the commutation 
relations = ±Rz and [R^,Rz] = ± 2 i?T^ Again, 

one performs an unitary transformation to the system’s 
Hamiltonian, U{t) = g,^i^R^+‘^phb^b)t^ represents it as 


where the Liouville superoperators Cqd and Cph describe, 
respectively, the QD’s and phonons’ dissipative effects. 
In the bare state representation QD’s dissipation pro¬ 
cesses are expressed by the following spontaneous emis¬ 
sion term: CqdP = - 7 [S'+,S'“p] - 'yc[Sz,Szp] + H.c. 

1^ . In the dressed-state basis and within the sec¬ 
ular approximation, i.e. 2D 3> 7 , the same processes 
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are described by three terms determined by the QD’s 
dressed-state decay rates: 7 + = 7 cost's + | 7 cSin^ ( 20 ), 
7 _ = 7 sin ‘*0 + |; 7 cSin^( 20 ) and 70 = |[ 7 sin^( 20 ) + 
7 cCOS^( 20 )]. Therefore, 

^qdP = -"f+[R'^, R~p] - "f-[R~, R'^'p] 

-lo[Rz,RzP] (8) 

The phonons from the multilayered acoustical cavity are 
allowed to interact with the environmental thermal reser¬ 
voir. In the rotating wave approximation this process is 
described by two terms corresponding respectively to the 
cavity’s damping and pumping effects at a rate propor¬ 
tional to K = uiphlQ determined by the cavity’s quality 
factor Q [s^: 

^phP = -k{1 + n)[b\bp] — Kn[b,b^ p] +II.C.. ( 9 ) 

Here, h is the mean thermal phonon number correspond¬ 
ing to frequency ujph and environmental temperature T . 

Once the master equation is determined by Eqs. dSD- 
m, it is solved by projecting the density operator firstly 
in the QD’s basis and then in the phonon field’s basis 
[s^. The projection in the QD’s dressed state basis leads 
after some rearrangements to a system of six coupled 
differential equations involving the following variables: 

P^^^ = P+++P—, P^"^^ = P++ - P—, 

p^^'> = - P-+b, = b'‘p+- + p-+b, 

= p+-b'^-bp_+, p^^'> = p+_b'^ + bp_+, ( 10 ) 

where p^J = {i\p\j), {i,j G |+), |-)}, are the QD’s 
density matrix elements. Finally, the projection in the 
phonon Fock states basis {|n)} gives a set of infinite dif¬ 
ferential equations, namely, 

r<'> = (py - pp) 

- 2K{l + n) - {n + l)pj]2i) 

- 2Kn ((n -b , ( 11 ) 

= _,5!^(p(3)+p(5)) 

- 2(7+-7-)pW-2(7++7_)Pi2) 

- 2^{1 + n) -{n + l)Pi%) 

- 2Kn[{n + l)Pi^^ -nPi%), ( 12 ) 

pm = (;,(.) _ pm _ pm _ p_m-j 

- I (/3(2n - 1) - (5) Pi^^ - ( 7 + + 7- + 47o)Pi'^ 

- ^{1 + n) (( 2 n - 1)PP - 2{n + l)P®i + 2 P,( 5 )^ 

- «:h((2n+l)Pf)_2nP,li\), (13) 


PW = _,(;3(2n-l)-(5)Pi3)-(^++^_+47o)pW 

- k{1 + n) ((2n - l)pW - 2(n + QP^^V + 

- «h((2n+l)pW_2^pD_)^^^ 

pO) = ('pa) ^ pa) _ pW^ ^ pO^)^^ 

- i (/3(2n -b 1) - (5) Pi®^ - {j+ -b 7 - + 47o)-Pi®^ 

- k{ 1 + n) ((2n + l)pW - 2(n + l)Pi%) 

- nn ({2n + 3)Pi®) - 2nPl^l, - 2Pi^^) , (15) 

Pi®) = -*(/3(2n+l)-5)Pi®)-(7++7-+47o)Pi®) 

- «(1 + n) ((2n + l)Pi®) - 2(n + l)Pi®+\) 

- nn ((2n + 3)Pi®) - 2nPi% - 2Pi^^) . (16) 

Here Pn^ = (n|p)*)|n) while S = ujph — 20 -b 2A. 

In the next Section we shall describe the cavity phonon 
dynamics in the steady-state via second-order phonon- 
phonon correlation function as well as the mean phonon 
number. 


III. RESULTS AND DISCUSSION 

Consequently, the average phonon number in the cav¬ 
ity mode is expressed as: 


(n) = (fetfe) = ^nPi^). (17) 

n—Q 


The nanocavity second-order phonon-phonon correlation 
function is dehned as usual [35|, 




{b^^bb) _ n(n - l)Pi^) 

( 6 t 6)2 („)2 


The system of Eqs. (HU) - (HU) as well as the infinite 
series in the expressions (H3 - HI must be truncated at 
a particular value of n = Nmax such that the variables of 
interest remain unchanged if one further increases N^ax 


In what follows, we shall study the system in the 
steady-state regime, i.e., Pn'^ = 0 for {z = 1 • • • 6}. The 
second-order correlation function given by Eq. HS and 
the average phonon number in Eq. HID are used to de¬ 
scribe the phonon field behaviors in the acoustical cavity 
mode [ 3 ^. Once truncated, the system of coupled equa¬ 
tions (HU) - HU) is solved by setting the model’s param¬ 
eters { 7 , 7 c, g, WpZt, K, h, 0} and A, respectively. 

A general overview over the steady-state system’s be¬ 
havior in the strong coupling regime or within the secu¬ 
lar approximation, i.e. when {A,/3 = 0}, is presented 
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FIG. 2: (color online) (a) The second-order phonon-phonon 
correlation function ( 0 ) (blue lines) and the mean phonon 
number in the cavity (n) (red curves) as a function of the QD- 
laser detuning A normalized by the Rabi frequency 2Q and be¬ 
yond the secular approximation (continuous lines), and within 
the secular approximation (dashed curves). Here n = 0.04, 
2 fI /7 = 25 and K /7 = 5 x 10“®. (b) (ji*-^^(0) as a function of 
the normalized damping rate K/'y for thermal baths at differ¬ 
ent temperatures and for 212/7 = 25, and A/{2fl) = —0.7. 
Here, from top to down one has n = 0.64, n = 0.16, n = 0.08, 
n = 0.04 and n = 0.01, respectively, (c) ^*-^^(0) (3D surface) 
and (n) (density plot) as functions of K/'y and 252/7. Here 
n = 0.04 whereas A/(252) = —0.7. The plot regions corre¬ 
sponding to > 1 and g^^\0) < 1 are represented in dif¬ 

ferent mesh styles. The plot region corresponding to (n) < 1 
is represented in gray color. The purple curves in hgures (b) 
and (c) are identical. Other parameters are; 7 c /7 = 0.1, 
fl /7 = 15, Ljph/^ = 35. 

in Fig. O The phonon field’s statistics is described 
by the second-order correlation function ^^^^( 0 ) with 
5 *-^^ (0) = 1 describing the Poissonian phonon distribution 
while when g^^^(O) < 1 one has sub-Poissonian phonon 
statistics. One observes that moderate laser-QD coupling 
strengths, i.e. 52, as well as lower temperatures give a 
more prominent sub-Poissonian phonon statistics. Fur¬ 
thermore, beyond the secular approximation, the phonon 
statistics exhibits quantum features, i.e. < 1 

(compare the corresponding curves in Fig. [SJi) and the 
effect is more pronounced for stronger QD-phonon cou¬ 
pling strengths. Fig. [^b) shows, beyond the secular 
approximation, the second-order phonon-phonon corre¬ 
lation function as a function of k/^. Here, again, around 
K > 7 one has a quantum phonon effect, i.e. sub- 
Poissonian phonon statistics. Thus, the contribution to 
the system’s dynamics of the fast rotating terms 
evaluated by the Eq. (O are essential for a sub-Poissonian 
quantum feature, see Fig.[51 In Fig. [3] we compare our re¬ 
sult within and beyond the secular approximation. The 


second-order correlation function estimated in both cases 
converges for higher and smaller cavity damping rates. 
However, for lower values of K /7 quantum features are 
proper only beyond the secular approximation (see the 
inset picture in Fig. [3]). Furthermore, the mean phonon 
number decreases in this particular case although it is 
quite high comparing to the k > 7 situation. 

IV. SUMMARY 

In summary, we have investigated the phonon quantum 
statistics in an acoustical nanocavity. A laser-pumped 
two-level quantum dot is embedded in the cavity con¬ 
tributing to phonon quantum dynamics. We have demon¬ 
strated that stronger QD-phonon-cavity coupling regimes 
lead to quantum features of the cavity phonon field in the 
steady-state. This QD-cavity interaction regime obliges 
to go beyond the secular approximation. Ignoring this 
fact would lead to an erroneous estimation of the phonon 
statistics for some parameter domains. 

Acknowledgement 

We are grateful to financial support via the research 
Grant No. 13.820.05.07/GF. 



Kh 

FIG. 3: (color online) The second-order phonon-phonon cor¬ 
relation function 5 f^^^( 0 ) (blue lines) and the mean phonon 
number in the cavity mode (n) (red curves) as a function 
of the cavity damping rate k normalized by the spontaneous 
emission rate 7 . The continuous curves are beyond the sec¬ 
ular approximation whereas the dashed ones are within the 
secular approximation. Here, 252/7 = 25, A/(2f2) = —0.7, 
n = 0.04, 7 c /7 = 0.1, g/'y — 15 and lOph/'y ~ 35. The con¬ 
tinuous blue line is identical with the purple lines of Fig. ©• 
The inset picture represents a close look on the behaviors of 
the second-order correlation functions in the regions around 
10"® < k/7 < 10"^ 
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